PASJ: Publ. Astron. Soc. Japan , -| (2008) 



The Light-cone Effect on the Clustering Statistics 
in the Cosmological Redshift Space 

Yasushi Suto, 1 ' 2 Hiromitsu Magira 1 

1 Department of Physics, University of Tokyo, Tokyo 113-0033 
E-mail(YS): suto@phys.s.u-tokyo. ac.jp 
2 Research Center for the Early Universe, School of Science, University of Tokyo, Tokyo 113-0033 
■ and 

Kazuhiro Yamamoto 

• Department of Physics, Hiroshima University, Higashi- Hiroshima 739-8526 

D ' 

r\ . (Received 1999 October 26; accepted 1999 December 14) 

Abstract 

We present a theoretical formalism to predict the two-point clustering statistics (the power spectrum 
and the two-point correlation function), simultaneously taking account of the linear velocity distortion, 
the nonlinear velocity distortion (hnger-of-god) , the cosmological redshift-space distortion and the light- 
cone effect. To demonstrate the importance of these effects in exploring the clustering of objects at high 
redshifts, we show several model predictions for magnitude-limited surveys of galaxies and quasars. This 
methodology provides a quantitative tool to confront theoretical models against the upcoming precision 
data on clustering in the universe. 
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Observations of clustering of the high-redshift universe suffer from a variety of contaminations including the linear 
redshift-space (velocity) distortion (Kaiser 1987; Hamilton 1998), the nonlinear redshift-space (velocity) distortion 
^ ■ (Davis & Peebles 1983; Suto & Suginohara 1991; Cole, Fisher, & Weinberg 1994, 1995), the cosmological redshift- 
• • , space (geometrical) distortion (Alcock & Paczyhski 1979; Ballinger, Peacock, & Heavens 1996; Matsubara & Suto 
' 1996), and the cosmological light-cone effect (Matsubara, Suto, & Szapudi 1997; Mataresse et al. 1997; Nakamura, 
K> ! Matsubara, & Suto 1998; Moscardini et al. 1998; de Laix & Starkman 1998; Yamamoto & Suto 1999; Suto et al. 
■ 1999; Nishioka & Yamamoto 1999). 

These effects have been already discussed in the previous papers quoted above, although separately and not in a 
unified fashion. In particular, Yamamoto, Nishioka, & Suto (1999) derived a formula for the power spectrum which 
takes account of the three effects except the cosmological redshift-space distortion, while Magira, Jing, & Suto (2000) 
extensively discussed the first three effects albeit neglecting the cosmological light-cone effect. The former approach 
requires to specify the (correct) values of both the density parameter £Iq and the cosmological constant Ao (see also 
Nishioka & Yamamoto 2000), and the latter analysis is applicable only for data over a narrow range of redshift (i.e., 
Az/z« 1). 

In this Paper, we develop for the first time the relevant formulae for two-point clustering statistics (the power 
spectrum and its Fourier transform, the two-point correlation function) which properly incorporate all the four 
effects mentioned above. To demonstrate the importance of these effects in exploring the clustering of objects at 
high redshifts, we present the predictions for future redshift surveys of galaxies and quasars, assuming spatially-flat 
cold dark matter (CDM) cosmology as a representative model. 



2. Formulae for the Two-point Clustering Statistics on a Light-cone in the Cosmological Redshift 
Space 

Our method to predict the two-point clustering statistics of objects on a light-cone in the cosmological redshift 
space proceeds as follows. First we have to specify the model of bias which relates the density field of objects, 
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<S(x, z), to that of mass <5 maS s(x,z). Our formulation presented here assumes the linear bias model (§2.1). While the 
realistic bias model is not well understood yet, the linear bias model is a fairly good approximation on large scales 
independently of the details of the biasing mechanism (e.g., Coles 1993, Matsubara 1999). Second, we take into 
account the redshift-spacc distortion induced by the peculiar velocity field (§2.2), using both linear theory (Kaiser 
1987) and a non- linear modeling for the finger-of-god effect (Peacock & Dodd 1996). Third, we correct for the 
anisotropy due to the geometry of the universe (§2.3). The resulting anisotropic power spectra are expanded in 
multipoles and translated to the corresponding moments of the two-point correlation functions (§2.4). Finally those 
multipoles defined at a given redshift are integrated over the past light-cone with the appropriate weight (§2.5), 
which yield the final expressions for the two-point clustering statistics of objects on a light-cone in the cosmological 
redshift space. 



2.1. model of bias 

Since we are mainly interested in the relation between the two-point statistics on a constant-time hypersurface in 
the real space and that on a light-cone hypersurface in the cosmological redshift space, we do not adopt a realistic, 
and consequently too complicated, model of bias. Rather we consider the case of the deterministic, linear and 
scale- independent bias: 

5(x, z) = b(z) <5 mass (x, z), (1) 

although the real bias might be scale-dependent, nonlinear and stochastic to some extent (e.g., Taruya, Koyama, & 
Soda 1999; Dekel & Lahav 1999; Matsubara 1999). Throughout the Paper, we explicitly use the subscript, mass, to 
indicate the quantities related to the mass density field, and those without the subscript correspond to the objects 
satisfying equation ([!]). 



2.2. redshift- space distortion due to the peculiar velocity 

The power spectrum distorted by the peculiar velocity field (neglecting the cosmological distortion for the moment), 
P( s >(k; z), is known to be well approximated by the following expression (Cole et al. 1995; Peacock & Dodds 1996): 



P (s) (k ± ,kfz) = b 2 (z)Pj^l{k;z) x 



-i 2 



D [fc||Op(aO] , 



(2) 



where k± and fcii are the comoving wavenumber perpendicular and parallel to the line-of-sight of an observer, and 
P^) s (k;z) is the mass power spectrum in real space. The second factor in the right-hand-side of equation (||) 
represents the linear redshift-space distortion derived by Kaiser (1987) adopting the distant-observer approximation 
and the scale-independent deterministic linear bias. The parameter f3(z) is defined by 



(3(z) 



1 dlnD + (z) 
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where D + (z) is the gravitational growth rate of the linear density fluctuations, and a is the cosmic scale factor. In 
the above, we introduce the density parameter, the cosmological constant, and the Hubble parameter at redshift z, 
which are related to their present values respectively as 



Q(z) 



Hp 

H(z) 



(i + z) 3 n , \(z) = 



H 



[H(z) 



Ao, 



H(z) = H ^n (l + zf + (l - fi„ - A )(l + zf + A . 
For definiteness we adopt an exponential distribution function for the pairwise peculiar velocity: 



fvivxi) 



1 



■ exp 
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which leads to the damping function in fc-space: 
1 



Z?[fc/i(7p] = 



l + (fc/!CT P ) 2 /2' 

with er P being the 1-dimensional pair-wise peculiar velocity dispersion. 



(4) 
(5) 

(G) 
(7) 
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2.3. redshift- space distortion due to the geometry of the universe 

Due to a general relativistic effect through the geometry of the universe, the observable separations perpendicular 
and parallel to the line-of-sight direction, x s± — (c/H )z58 and x s]] = (c/H )5z, are mapped differently to the 
corresponding comoving separations in real space x ± and xf 

x s± (z) = x ± cz/[H (l + z)d A (z)} = x ± /c ± (z), x H {z) = x\\H(z)/H = X\\/c\\[z), (8) 

with d A (z) being the angular diameter distance (Alcock & Paczyhski 1979; Matsubara & Suto 1996). The difference 
between c±(z) and c\\(z) generates the apparent anisotropy in the clustering statistics which should be isotropic in 
the comoving space. Specifically, the power spectrum in the cosmologicai redshift space, p( CRD ) 5 ig related to P' s ' 
defined in the comoving redshift space (eq.pl) as 



P^ CRD \k sU k sf z) = 



1 



c ± (z) 2 c\\{z) 



p(S) 



(9) 



c±{z) c||(z) : 

where the first factor comes from the Jacobian of the volume element dk 2 ± dk sfl , and k s ± = c±(z)k± and k s \\ — C|| (z)k\\ 
are the wavenumber (in the cosmologicai redshift space) perpendicular and parallel to the line-of-sight direction. 
Substituting equation (||), equation (^) reduces to 



P^{k s ^ k -z) 



b 2 (z) 



c ± {z) 2 c\\{z) 

1 -2 



_p(R) 
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~ 1 



2c 2 (z) 



(10) 
(11) 



where we introduce 

fc s = k 2 ± + fc^| , Hk = k s \\/k s , ?y = C||/c_L, 

(Ballinger, Peacock, & Heavens 1996; Magira, Jing, & Suto 2000). 

The two-point correlation function in the cosmologicai redshift space, (x s ±, x s f, z), is computed using equa- 

tion (10) as 



£ (crd) (x s ;z) = ^|p( CRD )(k s ;z)expHk s .x s )d 3 fc s 

; z) exp(— ik • x)d 3 fc 

lit) 6 J 



(27 

= ^ S \cxx s ±,c\\x s \\;z), 
where ^ s '(xx,xu; z) is the redshift-space correlation function defined through equation 



(12) 



2.4- multipole expansion 

Following Hamilton (1992), we decompose the power spectrum into harmonics: 



P(k,^ k ;z) = Pl(k)Li(fi k ), Pi{k;z) 



21 + 1 



l:even 



dfi k P{k,Hk;z)Li(ii k ), 



(13) 



where /i k is the direction cosine between the wavevector and the line-of-sight, and Li(/i k ) are the l-th order Legendre 
polynomials. Similarly, the two-point correlation function is decomposed as 



21 + 1 



lieven 



d/j, x £,(x,n x ; z)Li(n x ), 



using the direction cosine fi x between the separation vector and the line-of-sight. 
The above multipole moments satisfy the following relations: 

[°° 1 
Pi(k; z) = Am 1 I &(x; z)ji(kx)x 2 dx, £i(x; z) 



2n 2 i l 



Pi(k; z)ji(kx)k dk, 



(14) 



(15) 



with ji(kx) being the spherical Bessel functions. Substituting P^ CRD \k s , ^ k ; z) and ^ CRD ^(x s ; z) in equations (|T 
to (p^|), one can compute the moments, Pi CRD \k s ; z) and ^l° RD \x s ; z) at a given z. Actually the above relations 
are essential in predicting the two-point correlation function in the cosmologicai redshift space with the nonlinear 
effects. 
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2.5. the light-cone effect 

Finally we incorporate the light-cone effect following the expression by Yamamoto & Suto (1999). Since 
P[ C (k s ', z) and 6 (A; s ; z) are defined in redshift space, the proper weight should be 

d 3 s (cR D ) ^ (z)n CR D(z)] 2 = d 3 x ^ z ) n ^^2 c± ( z )2 C||(z)j ( 16 ) 

where n§ RD (z) and nQ om (z) denote number densities of the objects in cosmological redshift space and comoving space, 
respectively, and 4>{z) is the selection function determined by the observational target selection and the luminosity 
function of the objects. Then the final expressions reduce to 

dz^ [4>(z)n—{z)fc^zfc^z)P^°\k s] z) 

p ( L C,CR D)(M = J Znm _g (17) 

dz^- [<j){z)nir(z)] 2 ct{zfc^z) 

dz d lL [^z)n—{z)fc^zfc^z)^(x s -z) 
#°' CRD) (^) = fZ l Z rlv , (18) 



dz^r [(f){z)n'- o ' n {z)] 2 c^{z) 2 c^z) 



where z max and z m i n denote the redshift range of the survey, dV c /dz = d 2 ,(z) / H \z) is the comoving volume element 
per unit solid angle. 

Note that k s and x s defined in P^ CRO \k s ; z) and ^ CRD (a; s ; z) are related to their comoving counterparts at z through 
equations ( fi~ij ) and (||) while those in pj ljC '' CRD \k s ) and £,l LC ' CRD \x s ) are not specifically related to any comoving 
wavenumber and separation. Rather they correspond to the quantities averaged over the range of z satisfying the 
observable conditions of x s = (c/Hq)\/Sz 2 + z 2 59 2 and k s — 2ir/x s . 



3. Predictions for Galaxy and QSO Samples in CDM Models 

As in Yamamoto, Nishioka, & Suto (1999), we consider SCDM (standard cold dark matter) and LCDM (Lambda 
cold dark matter) models, which have (fio, Ao, ft, eg) = (1-0,0.0,0.5,0.6) and (0.3,0.7,0.7,1.0), respectively. These 
sets of cosmological parameters are chosen so as to reproduce the observed cluster abundance (Kitayama & Suto 
1997). We use the fitting formulae of Peacock & Dodds (1996) and Mo, Jing, & Borner (1997) for the nonlinear 
power spectrum Pj^} s (k] z) and the peculiar velocity dispersions <7 P , respectively (see also Suto et al. 1999; Magira 
et al. 2000). Also we take into account the selection functions relevant to the upcoming SDSS spectroscopic samples 
of galaxies and quasars adopting the B-band limiting magnitudes of 19 and 20, respectively. Unlike Yamamoto & 
Suto (1999) and Yamamoto, Nishioka, & Suto (1999), we properly transform the QSO luminosity function originally 
fitted assuming the Einstein - de Sitter universe (Wallington & Narayan 1993; Nakamura & Suto 1997) to that in a 
given cosmological model. 

Figure 1 compares the predictions for the angle-averaged (monopole) power spectra under various approximations. 
Upper and lower panels adopt the selection functions appropriate for galaxies in < z < z max = 0.2 and QSOs 
in < z < z m ax = 5, respectively. Left and right panels present the results in SCDM and LCDM models. For 
simplicity we adopt a scale-independent linear bias model of Fry (1996): 

b(z) = l + -^ T -[b(k,z = 0)-l], (19) 

with b(k, z = 0) = 1 and 1.5 for galaxies and quasars, respectively. We present the results normalized by the real- 
space power spectrum in linear theory (Bardeen et al. 1986); P{f \k; z — 0), P^ S \k; z = z max ), PQ CRD \k s ; z = z max ), 
and PQ LC ' CRD ^(fc s ) are computed using the nonlinear power spectrum of Peacock & Dodds (1996), and plotted in 
dashed, dash-dotted, dotted and solid lines, respectively. 

Consider the results for the galaxy sample (upper panels) first. On linear scales (k < O.lftMpc -1 ), Pq (fc;2 = 0) 
plotted in dashed lines is enhanced relative to that in real space mainly due to the linear redshift-space distortion 
(The Kaiser factor in eq.j2|). For nonlinear scales, the nonlinear gravitational evolution increases the power spectrum 
in real space while the finger-of-god effect suppresses that in redshift space. Thus the net result is sensitive to the 
shape and the amplitude of the fluctuation spectrum itself; in the LCDM model that we adopted, the nonlinear 
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Fig. 1.. The light-cone and cosmological redshift-space distortion effects on angle-averaged power spectra. Upper and lower panels 
correspond to magnitude-limited samples of galaxies (B < 19 in < z < z m ax = 0.2; no bias model) and QSOs (B < 20 in 
< z < ,z m ax = 5; Fry's linear bias model), respectively. Dashed: (k; z = 0)/P( n,1Jn )(fc;2 = 0), dash-dotted: Pp S '(fc;2: = 2 max ), 
dotted: P^ CRD \fc s ; z = z max ), solid: Pq L °' CRD ' 1 (fc a ). 



gravitational growth in real space is stronger than the suppression due to the finger-of-god effect. Thus P ' (fc; z = 0) 
becomes larger than its real-space counterpart in linear theory. In the SCDM mode, however, this is opposite and 
P (k;z = 0) becomes smaller. 

The power spectra at z = 0.2 (dash-dotted lines) are smaller than those at z = by the corresponding growth 
factor of the fluctuations, and one might expect that the amplitude of the power spectra on the light-cone (solid 
lines) would be in-between the two. While this is correct if we use the comoving wavenumber, the actual observation 
on the light-cone in the cosmological redshift space should be expressed in terms of A; s (eq.[(0|). If we plot the 
power spectra at z = 0.2 taking into account the geometrical distortion, Pq (As a ; z = 0.2) in dotted lines becomes 
significantly larger than P^ {k; z = 0.2). Therefore Pq C ' GHD ' (fc s ) should take a value between those of Pg GHD ^(fc s ; z = 
0) = P^\k] z = 0) and Po° RD ^(fc s ; z = 0.2). This explains the qualitatively features shown in upper panels of figure 
1. As a result, both the cosmological redshift-space distortion and the light-cone effect substantially change the 
predicted shape and amplitude of the power spectra, even for the galaxy sample (Nakamura, Matsubara, & Suto 
1998). 

The results for the QSO sample can be basically understood in a similar manner except that the evolution of bias 
makes significant difference since the sample extends to much higher redshifts. Although it is true that the results 
are sensitive to the model of bias, we believe that those on linear scales are fairly insensitive to the complexities of 
the bias; it has been shown that the bias on such scales can be well approximated by the linear and scale-independent 
model (Coles 1993; Matsubara 1999; Magira, Taruya, Jing & Suto 2000). 

Figure 2 shows the results for the angle-averaged (monopole) two-point correlation functions, exactly corresponding 
to those in figure 1. The results in this figure can be also understood by the analogy of those presented in figure 1 at 
k rsj 2t:/x. Unlike the power spectra, however, two-point correlation functions are not positive definite. The funny 
features in figure 2 on scales larger than 30/i _1 Mpc (100/i _1 Mpc) in SCDM (LCDM) originate from the fact that 
£(R.,im) ^ 2 = 0) becomes negative there. 

In the linear and deterministic bias model that we adopted here, the power spectra of objects at z in redshift space 
(without the geometrical effect) are proportional to the mass power spectrum in real space, 

p $£L(k;z)- Hamilton 
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Fig. 2.. The same as Figjljon angle- averaged two-point correlation functions. 



(1992) showed that in linear theory of redshift-space distortion (Kaiser 1987), the quadrupole to monopole ratio 
provides a measure of estimating /3 through the following relation: 

4„, , 4 



C| ,lm (»; z) 



(20) 



3 / £ Q ' m (xw; z)w 2 dw 



Taking account of the nonlinear redshift-space distortion, however, the ratio becomes scale-dependent; our present 
model (||) with equation (0) yields 



pS(k;z) 
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, B{n) = — 
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arctan(^) 



^ 3 3arctan(«;) 



(21) 
(22) 



with k(z) = ka P (z)/y2Ho (Magira, Jing, & Suto 2000; Yamamoto, Nishioka, & Suto 1999; see also Cole, Fisher, 
& Weinberg 1995). We note here that Pf(fc;z) is n °t an y m ore positive definite and in fact changes sign at 
cosmologically interesting scales (see figure 3 below). 

The above ratio of monopole and quadrupole turns out to be independent of the underlying mass power spectrum, 
which is regarded as one of the advantages of estimating the /3-parameter using the redshift-space distortion. The 
cosmological redshift-space distortion (eqs. @ and @) and the light -cone effect (eqs.(l7) and |Tq]), however, further 
complicate the behavior of the final predictions for the two-point statistics, which are plotted in figure 3. In fact, 
even the shape of Pq LC,crd ' 1 (fc s ) and p 2 ( LC,CRD ) (fc s ) f or QSOs is sensitive to the bias model and the cosmological 
parameters. 

For the galaxy sample, P 2 (LC ' CRD) (fc s )/P (LC,CRD) (£; s ) in large scales approaches 1.06 (0.65) for SCDM (LCDM), 
which should be compared with 50/49 (0.57) in linear theory (eq.^cj). This implies that such estimated values of 
(3 are 1.05 (0.57), in contrast to the true values of 1.0 (0.49). Thus even for the shallow samples of galaxies, the 
present effects systematically changes the estimate of (3 by (5 ~ 20)% level (Nakamura, Matsubara, & Suto 1998). 
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k a [hMpc" 1 ] x s [h-'Mpe] 

Fig. 3.. Predictions for the two-point statistics on a light-cone in the cosmological redshift space. The results in SCDM and LCDM 
models are plotted in thin and thick lines, respectively. Left: monopolc, quadrupole and quadrupole-to-monopole ratio of the 
power spectrum (from top to bottom). Right: monopole, quadrupole and quadrupole-to-monopole ratio of the two-point correlation 
functions (from top to bottom). 

4. Conclusions and Discussion 

We have developed for the first time a theoretical formalism to predict the two-point clustering statistics on a 
light-cone in the cosmological redshift space, p( LC > CRD ) (fc s ) and ^Q LC,CRD ' l (a; s ), which were only discussed separately 
before (Yamamoto, Nishioka, & Suto 1999; Magira, Jing, & Suto 2000). This methodology fully takes into account 
all the well-defined physical, i.e., gravitational and dynamical effects, and enables one to quantitatively confront the 
future precision data on the clustering statistics against the cosmological model predictions modulo the uncertainty 
of the bias. 

In fact, since the resulting predictions are sensitive to the bias, which is unlikely to quantitatively be specified by 
theory, the present methodology will find two completely different applications. For relatively shallower catalogues 
like galaxy samples, the evolution of bias is not supposed to be so strong. Thus one may estimate the cosmological 
parameters from the observed degree of the redshift distortion as has been conducted conventionally. Most impor- 
tantly we can correct for the systematics due to the light-cone and geometrical distortion effects which affect the 
estimate of the parameters by ~ 10% (§3). This is generally the case on linear and quasi-nonlinear scales where it 
has been shown that the bias on such scales can be well approximated by the linear and scale-independent model 
(Coles 1993; Matsubara 1999; Magira, Taruya, Jing & Suto 2000) as we have assumed in this paper. 

Alternatively, for deeper catalogues like high-redshift quasar samples, one can extract information on the nonlin- 
earity, scale-dependence and stochasticity of the object-dependent bias (Taruya, Koyama, & Soda 1999; Dekel & 
Lahav 1999; Blanton et al. 1999) only by correcting the observed data on the basis of our formulae. In this case one 
should adopt a set of cosmological parameters a priori, but they will be provided both from the low-redshift analysis 
described above and from the precision data of the cosmic microwave background and supernovae la. 

In a sense, the former approach uses the light-cone and geometrical distortion effects as real cosmological signals, 
while the latter regards them as inevitable, but physically removable, noises. In both cases, the present methodology 
is essential in properly interpreting the observations of the universe at high redshifts. 

We thank Shin Sasaki for suggesting the proper scaling of the QSO luminosity function to an arbitrary cosmological 
model. This research was supported in part by the Grants-in-Aid by the Ministry of Education, Science, Sports and 
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Culture of Japan to RESCEU (07CE2002) and to K.Y. (11640280), and by the Inamori Foundation. 
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